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Abstract 

Little Higgs models offer a new way to address the hierarchy problem, and give 
rise to a weakly-coupled Higgs sector. These theories predict the existence of new 
states which are necessary to cancel the quadratic divergences of the Standard Model. 
The simplest version of these models, the Littlest Higgs, is based on an SU{5)/SO{5) 
non-linear sigma model and predicts that four new gauge bosons, a weak isosinglet 
quark, F, with Q = 2/3, as well as an isotriplet scalar field exist at the TeV scale. We 
consider the contributions of these new states to precision electroweak observables, and 
examine their production at the Tevatron. We thoroughly explore the parameter space 
of this model and find that small regions are allowed by the precision data where the 
model parameters take on their natural values. These regions are, however, excluded 
by the Tevatron data. Combined, the direct and indirect effects of these new states 
constrain the ‘decay constant’ / ^ 3.5 TeV and mf ^ 7 TeV. These bounds imply 
that significant fine-tuning be present in order for this model to resolve the hierarchy 
problem. 
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1 Introduction 


The Standard Model (SM) of electroweak physics is a remarkable achievement. Precision 
electroweak experiments have probed the SM at the level of qnantum effects and have con¬ 
firmed every featnre of the theory. In particular, the set of precision electroweak (EW) data 
has tested the SM beyond one-loop level and no significant deviations from SM predictions 
have been observed. The symmetry breaking sector of the SM has been investigated through 
its virtual effects and precision measurements strongly prefer the existence of a weakly cou¬ 
pled Higgs boson by constraining its mass to be mriH < 193 GeV at 95% CL jT]. This upper 
bound may be relaxed if certain classes of new physics contributions beyond the SM are 
present at the TeV scale and compensate for the indirect effects of a heavier scalar sector 
[2j. However, because of decoupling, models with an elementary Higgs scalar typically yield 

small contributions to the electroweak radiative corrections. Direct evidence for the elec¬ 
troweak symmetry breaking dynamics has yet to be observed, and searches for Higgs boson 

production have placed the compatible lower bound rriH > 114.4 GeV. 

The existence of a weakly coupled Higgs sector generates the hierarchy problem, 
which is a long-standing puzzle in particle physics. Few candidates exist for the mecha¬ 
nism which solves this problem. Supersymmetry naturally stabilizes the hierarchy, as the 
quadratically divergent loop contributions to the Higgs mass cancel between the fermion and 
boson contributions, provided the supersymmetry breaking scale is near a TeV [S]. Novel 
theories with extra dimensions exploit the geometry of the higher dimensional spacetime to 
resolve the hierarchy |3] and result in numerous phenomenological consequences at the TeV 
scale jH] . Supersymmetry has the additional feature of remaining weakly coupled up to high 
scales, whereas quantum gravity becomes strong at the TeV scale in these extra-dimensional 
models. Experimental confirmation of either scenario has yet to occur. 

Recently, another scenario has been developed which attacks the hierarchy problem 
in a new way, while maintaining a weakly coupled scalar sector. This scenario is known 
as the Little Higgs inmi where the Higgs is effectively a pseudo-Goldstone boson. These 
theories are realizations of earlier attempts [HI to stablize a Higgs which arises as a pseudo- 
Goldstone boson resulting from a spontaneously broken non-linear approximate symmetry. 
This approximate global symmetry protects the Higgs vev relative to the UV cut-off of the 
theory. A, which appears at a higher scale. The recent progress was attained in models of 
dimensional deconstruction |H], where the quadratically divergent corrections to the Higgs 
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mass were shown to cancel with contribntions from new helds in non-trivial descriptions of 
theory space. The phenomenology of these models was examined in HU- The recent Little 
Higgs models have the advantage of not reqniring a non-trivial theory space in order for the 
qnadratic divergences to be eliminated. The novel feature of these cancellations, is that the 
divergent contributions from a particular particle are cancelled by a new particle of the same 
spin. Hence, these theories predict the existence of new Q = -1-2/3 quarks, gauge bosons, 
and scalars, all with masses at the TeV scale, in order to remove the relevant divergences. 
In addition, a small mass for the Higgs boson, of order 100 GeV, is naturally obtained from 
multi-loop corrections. 

The most economical model of this type to date, known as the Littlest Higgs 0, is 
based on an SU{5)/SO{5) nonlinear sigma model. The UV cut-off of this theory occurs 
at A dvr/ 10 TeV, where the nonlinear sigma model becomes strongly coupled. Here 
/ is the decay constant of the pseudo-Goldstone boson and is necessarily of order a TeV. 
The SU{5) SO{5) symmetry breaking is generated at the scale A by strong interactions 

similar to technicolor which act only at the high scale A and have no remnants at a TeV. The 
14 Goldstone bosons remaining after this symmetry breaking yield a physical doublet and 
a complex triplet under SU{2)] the remaining helds are eaten by a Higgs-like mechanism 
when the intermediate symmetry is broken. The one-loop contributions to the quadratic 
divergences from the electroweak gauge sector are removed by making use of the SU{5) 
weakly gauged subgroup [SU{2) x f/(l)]^. The divergent contributions must then involve 
couplings of both groups and hence hrst appear at two-loop order. This combination of 
gauge couplings breaks the global symmetries and generates a small mass of order 100 GeV 
at the two-loop level for the scalar doublet held. At one-loop order, the Higgs mass is not 
sensitive to the ehective high scale A, and its mass is thus relatively stable assuming A ~ 10 
TeV. The scalar triplet acquires a mass at the TeV scale from one-loop gauge interactions. A 
massive vector-like SU{2) singlet Q = -1-2/3 quark is responsible for canceling the one-loop 
divergences originating from the SM top-quark. 

The physical spectrum of this model below a TeV is thus simply that of the SM with 
a single light Higgs. At the TeV scale four new gauge bosons V (an electroweak triplet and 
singlet) appear, as well as the scalar triplet cp, and a single vector-like quark t'. For the 
quadratically divergent corrections to cancel naturally, and not be hne-tuned, the mass of 
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the new quark is constrained to be [H] 


mti ^ 2 TeV 


mn 

.200 GeV. 


which in turn implies 


/ < 1 TeV 


ruH 


1 2 


200 GeV 


In addition, naturalness also implies 


mvi < 6 TeV 


rriH 

.200 GeV. 


7710 < 10 TeV. 


( 1 ) 

( 2 ) 

( 3 ) 


The Littlest Higgs, based on a SU{5)/SO{5) non-linear sigma model, is the simplest 
scenario of this type in the sense that it introduces a minimal number of new helds. Other 
models have been investigated 0 and are based on the cosets SU{k)/Sp{k), SU{kY/SU{k), 
and SU{kY/SU{k)"^, where k is the number of strongly interacting fermions at the scale A. 
These cases differ from the Littlest Higgs in that they require the introduction of additional 
scalars, gauge bosons, and vector-like quarks at the TeV scale, and sometimes extra light 
scalar helds can arise at ~ 100 GeV. 

In this paper, we examine the effects of the components of the Little Higgs on precision 
electroweak observables and in direct searches for new particles. We concentrate on the 
minimal model here, the Littlest Higgs [H], as it contains all the essential ingredients and 
provides a representative example of these types of theories. We would expect the more 
complicated scenarios listed above to be even further constrained by experiment in the 
absence of any hne-tuning due to the presence of the additional new particles at the TeV 
scale. Gonsistency with the global precision electroweak data set is a tough barrier for any 
new model to pass, and we hnd that the Littlest Higgs does not necessarily decouple quickly 
enough for most of the parameter space. A thorough examination of the parameter space 
reveals small regions where precision data, alone, allows for scales in the theory to take on 
their natural values, i.e., f ^ 2 TeV. However, we hnd that the most stringent constraint 
arises from the direct search for new gauge bosons in Drell-Yan production at the Tevatron 
CH and that this excludes these regions. Together, the direct and indirect data samples place 
the constraints / ^ 3.5 — 4.0 TeV, and ^ 7 TeV, indicating that signihcant hne-tuning 
is required if this scenario is to work. 
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The outline for the remainder of this paper is as follows: in the next section we 
present the formalism for our analysis. We start with a description of the parameters in 
this theory and then give a derivation of the shifts that occur in the precision electroweak 
observables. In section 3, we present our numerical results. We hrst examine the bounds 
obtained from direct searches for new gauge bosons at the Tevatron. We then study the case 
where the new particles take on SM values for their coupling strengths and perform a £t 
to the electroweak data set and obtain tight constraints on the model. We then vary these 
couplings within natural ranges and demonstrate that there is a region of parameter space 
where the precision data bounds are relaxed, but is still excluded by the Tevatron data. The 
hnal section contains our conclusions. 


2 Formalism 


To begin our analysis, we follow the notation of Arkani-Hamed et al.^ and generate the 
masses of the [SU{2) x gauge bosons through the leading two-derivative kinetic term 

for the non-linear sigma model 

£ = (4) 

where / ~ 1 TeV and the S = exp{2ill/f)Tio represents the Goldstone boson helds. The 
fourteen helds 11 include the SM weak iso-doublet h, an isotriplet 0 as well as the true 
Goldstone bosons which are eaten in the breaking of the intermediate symmetry to the SM. 
We choose the normalization 


n 


^ ht /2 01/^2 ^ 

V2 h*/2 

^ 0 / n /2 h ^/2 ^ 


( 5 ) 


such that the normalization of the kinetic term above is in agreement with that in earlier work 
on the non-linear sigma model ^2]- fa addition to the usual Higgs doublet vev, the triplet 
(j) also acquires an SU{2) breaking vev. We demonstrate in the Appendix that including the 
contribution from the triplet vev to the precision electroweak observables only strengthens 
our conclusions; here we adopt a conservative approach and neglect this piece in the following. 
We note that the presence of a scalar triplet is not a generic feature of Little Higgs models. 
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The covariant derivative of S is given by 


DE = as - iQkWmp: + SQf) - ig'MYj: + SF,^), (6) 

where a sum over the index /c = 1,2 is understood and the ‘charge’ (Q“) and ‘hypercharge’ 
(F) matrices are given in Ref. jB]. To proceed further in the calculation of the gauge boson 
masses it is useful to define the following combinations of the above [SU{2) x U{l)]k gauge 
couplings: 


9 


9t 


9192 


91 + 92 


9l 


9i + 92 


9 = 


9t = 


9192 


9? + 92 


9{ 


^9i +92 


Inverting these relations yields 


91 = (9^ + g?)‘'^ 

92 = ^(/ + 9?)‘-'^ 

9t 


( 7 ) 


( 8 ) 


and similarly for g[ 2 - To zeroth-order in the Higgs vev, one linear combination of the gauge 
bosons obtains a mass at the TeV scale when the intermediate symmetry is broken, while 
the second set remains massless. The massive neutral fields and their corresponding masses 
are given by 


Bh 

Zh 


g[Bi - g'^B^ 
^9?+ 92 



9iWl - g2Wi 
\j9i+ 92 


m 


2 1 {9^ + 9lf 

2 gl 



(9) 


The corresponding charged isospin partner of the F/j, Wj^, obtains an identical mass: 




9iWt - 92Wt 


9i+92 




^{9^ + 9l;r 


9t 


f 


( 10 ) 
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The orthogonal linear combinations 


fi(0) = 

92^1 + 91 B 2 

I 5 


^9i +92 

pj/3(0) ^ 

g2Wf + giWl 

J 9 1 + 9 2 


( 11 ) 


together with the corresponding charged field, remain massless at this order. They 

will play the roles of the SM fields and acquire masses through the usual Higgs mechanism 
which are proportional to the Higgs vev, u. 

In performing the comparison to precision electroweak and collider data, one must 
include the shifts in the SM-like gauge boson masses due to mixing with the heavier states, 
i.e., the Higgs vev induces a finite mixing between the states Bh, Zh, Wh and the correspond¬ 
ing lighter fields. For the neutral fields, which form a 4 x 4 mass matrix, this mixing is 
most easily examined via a power series expansion in the ratio of the Higgs vev relative to 
/. For the charged fields, since the mass matrices are now only 2x2, exact eigenvalues in 
a compact form can be easily found. The leading corrections to the masses of the lighter 
SM-like fields are of relative order 5^ = (z^//)^ ~ 0.01 so that, numerically, corrections 
can be safely neglected. When discussing the effects of the heavy gauge fields on precision 
measurements, we will consistently keep all terms of relative order 5^. Turning now to the 
light mass eigenstates, we define the zeroth-order weak mixing angle as follows: 



so that the physical 7 and Z bosons become (to O (5^)) 

H = ^ 

Z = c®W3(o) ^ ^(0)^(0) ^ ^ ^ (^3) 

where Sb^ , defined below in Eq. 15. For the physical bosons we obtain 




\/9? + I 4(9? + 9l) / ym 



g| (g? - al) ,2! 

4(9? +sit J 


Wo" 


(14) 
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which can easily be written in terms of the unmixed fields and 

_ 9i92{gi - 92 ) ^2 ^l± QgN 

4(^2 + ^ 2)2 I 

Note that the mixing between these states vanishes in the limit gi = g 2 or, more specifically, 
when gt = g- The physical Z contains a small admixture of the heavy fields proportional to 
^Bhi^Zh- These parameters, which define the mixing between the light SM-like Z and the 
heavy neutral gauge bosons, are given to order 5^ by: 


= 


b9t [91-9"^) \j9^ + 9'^ 
^ (9'^ + 9t^y 


4 (^2 + ^2)2 1 1 

Again we note that these mixing terms vanish in the symmetric limit when gt = g and 
g'l- = g'. From this we see that all light-heavy gauge boson mixing vanishes in the symmetric 
case. This will have important consequences below. The masses of the physical and Z 
bosons are then given by (to O (5^)): 

,,, 2(0) f, ,2 5 i9t + 9^) {g't^ -gS? + V 7^'^' \ 


Ml = 

2 (0) 

Mly = 

2 (0) 

where 



5 {gt + 9^) - gS^ + 9^) + ^9^91 [g^^ + gf) - iQg^^gh? 

Q {g"^ + 9tf {g'"^ + 9?)^ 


^2 bg^ - 2g^gl + bgt 
24 {g^ + gh" 


2(0) _ 1 / 2 I '2\ 2 2 ( 0 ) _ ^22 /lo 

- ^[9 + 9 ~ ^9 ^ 

are the usual gauge boson masses that satisfy the tree-level SM relations: . 


2 ( 0 ) _ 1 2„2 


^ = ^9 


These shifts in the masses of the gauge fields contribute to the deviation of the p- 
parameter from its SM value. These contributions alone yield 

_ 6M^ 5Mz . . 

^PG 2 ( 0 ) 2 ( 0 ) 


5 (9?-g'^ 


8 \9t +9 
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Note that as we might have anticipated the gauge contribution to 6p vanishes in the sym¬ 
metric limit, he., when g[ = g' . This is expected as in this limit the W and Z are just 

the usual SM helds. As we will discuss below there is an additional new contribution to 
5p beyond that arising from the tree level shifts of the gauge boson masses described here. 

This is due to the mixing of the top quark with the new vectorlike isosinglet field, i!. This 
contribution, which we denote by Sptop, arises from the new t and t' quark loops that do not 
appear in the SM. As we will see, these mixing induced contributions can be numerically 
signihcant in some regions of the parameter space and must be included. The reason for 
their signihcance is that such terms can be parametrically enhanced by powers of 
or (though correspondingly suppressed by powers of the t — t' mixing angle). Note 

that even though the additional isosinglet quark is vector-like, since it mixes with the top 
quark it does not fully decouple, he., it can lead to small but important corrections to 5p via 
mixing over some parameter space regions-particularly those where the t — t' mixing angle 
is large. (Here we do not mean that as the t' mass increases its contribution to 5p does not 
go to zero, he., the classical meaning of decoupling.) Thus the complete new contribution 
to the shift in the p-parameter is the sum of Sptop and Spa- 

To proceed further with an analysis of precision data and collider constraints, we must 
examine how the SM fermions transform under the group Gi x G 2 with Gk = [SU{2) xU{l)]k- 
Here we take the fermions to transform non-trivially only under Gi. This choice is not unique 
and at least two other possibilities exist. The hrst alternative is when the fermions transform 
identically under both Gi and G 2 . This possibility is excluded here since it would require 
a doubling of the number of fermions. A second somewhat convoluted alternative would 
have the fermions transform under, say, Gi as well as the U{1) piece of G 2 . In this case, 
doubling of the number of fermions is not required, but such a possibility would allow for an 
additional free parameter in the fermion couplings to the Z. It would seem that our choice 
is more natural and economical and avoids the issue of fermion doubling. However, we will 
make some particular notes in our numerical analysis below of the dependency of the choice 
of fermion couplings to the heavy hypercharge gauge boson. 

The t and t' quarks contribute at one-loop to 6 p as discussed above, as well as to 
Rb, the ratio of the b-quark to total hadronic width of the Z. In order to compute these 
contributions we must examine the t — t' mass matrix which can be obtained from the 
Yukawa couplings provided in j^]. The diagonalization of this matrix requires a bi-unitary 
transformation, which for real Yukawa couplings implies that two separate rotations on the t 


and t' fields, Ol,r, must be performed. Here Ol,r rotates the left(right)-handed components 
of the quark helds. It is important to note that the rotation Or is actually unphysical since 
both tji and transform identically under Gi. Ol, however, is physical and leads to, e.g., 
flavor violating Ztt' couplings. Since the mass matrix is a simple 2x2 rotation, the matrix 
Ol contains only a single parameter: a mixing angle 6t. Writing the ratio of the Yukawa 
couplings as r = X\/X\ and recalling 5 = u/f one hnds 


from which it follows that 

r = V 25 / tan(20i) - 1 + 5‘^/2 . (21) 

Since by dehnition r > 0, this immediately provides a bound on 9t. We hnd that 

tan(0i) < ^ . (22) 

Since 5 is relatively small this greatly restricts the range of Of. The eigenvalues of the t — t' 
mass matrix are now easily obtained with the smaller one being identihed with the physical 
top quark. Given 6,9t and the hxed experimental value of rrit = 174.3 ± 5.1 GeV [12], the 
mass of the t', rrit', is now completely determined: 


ml _ I+ r + 6^2 + ^(1 + r + 6^2)^ - 2rd^ 

1 + r + <572 - ^(1 + r + 572)2 - 2 r (52 ‘ 


(23) 


2 f2 


ml is then given by the numerator in the above equation times a factor of Xff 


For a given set of values for the parameters {S,9t), the contribution to the shift in 
the p-parameter from the fermion sector can now be fully calculated. (For simplicity we set 
Vtb = 1 in what follows.) The couplings of the t and t' to the W and Z helds are altered 
by the presence of the mixing and are represented by the parameter 9^, e.g., ordinarily 
the isosinglet t' would not couple to the W. We hnd Wib ~ q while Wi'b ~ St where 
Ct[st] = cos(6*t) [sin(6*t)]. Note that to this order in 6 we can neglect the ehects of gauge boson 
mixing on the couplings of the W and Z. For the Z, the right-handed couplings of the t and 
t' are not modihed but the left-handed couplings are now given by 


at = 


8G/M| 


■ci/2 


2x, 


-/3) 
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( 24 ) 



SG/Ml 

[sl/2 - 2x^/3) 

72 ' 

9^ = 9^ = ^ 






where = sin^ 9]v can be taken to be the on-shell valne to this order in 6. Note the presence 
of the flavor-changing conpling that we allnded to above. There are now five graphs which 
contribnte to the vacnnm polarization of the Z boson and two graphs for the W] the usnal 
SM contribntion must be subtracted from these new contributions to Sptop- To proceed let 
us write the general couplings of a pair of fermions to a gauge boson as 

^ + a'^b)f 2 X^ + h.c. (25) 

In this language, we now find that 

Sp =- SzCw), (26) 

where cw = cos 6w, which, again to this order, can be taken to be the on-shell SM value, 
and 

6x = 2[2mim2 {v‘^ — a^) — {mf + + a^)]\og{mim 2 /— 2mim2{v'^ — a^) (27) 

+ (^2 ~ rril)~^[ 2 mim 2 {rril + rnl){v‘^ — a^) — {m\ + rn 2 ){v‘^ + a^)] log(m 2 /m^), 

with mi ,2 being the relevant fermion masses. Here, p, is an arbitrary renormalization scale 
that cancels from the hnal expression after all the individual diagrams are summed over. 
For the Z boson vacuum polarization diagrams, the intermediate states which contribute 
are tt, i't', tt', t't and bb, while for the W they are ib and t'b, and their couplings are as given 
above. Using this expression we can now sum the new contributions and subtract those of 
the SM thus obtaining the contribution to the shift in the p-parameter from the fermion 
sector of this model. 

We now have the necessary ingredients to calculate the shifts in the precision elec- 
troweak observables. We perform our calculations in the on-shell renormalization scheme 
and begin by examining the effects related to the dehnition of Gf in this model. In the on- 
shell renormalization scheme, the boson mass is dehned through Gf, which is derived 
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from muon decay, and a, the fine structure constant. Note that in the on-shell scheme, Gf 
is employed as an input parameter. Several effects now modify this relation: (1) there are 
contributions from the exchange of the new gauge boson in /i-decay; (2) the coupling 
of fermions to the physical are shifted from the expression in Eq. [7| due to mixing with 
the as presented in Eq. [T3| 

gphysicai = gw = g{^- ^g) , ^g = ’ ( 2 ^) 

4 (gi+gi) 

(3) the expression for g in terms of a becomes g^ = 47ra/s|y^, where is defined in Ea. 1121 
this must be written in terms of the physical and Z masses. After implementing these 
changes, we arrive at the following equation relating Mw to Gf and a 

=_ TO (1 - 

^/2 2M^ (1 - (1c^(5p/s^) ' 

Here 6 p = Spc + Sptop- Solving this for to O (5^), we hnd the following shift in the 
mass from its predicted SM value. 


SM^ = 


1 ~ 




g^g^ 


^w 


2 ( 9 ^ + a'ff 



(30) 


We define {Mw)sm as the mass of the W boson given by the SM when the Higgs boson 
mass, rriH, as well as Mz, Gf, a{M‘^), as{M‘^) and rrit are used as input, and when all of the 
usual SM radiative corrections are incorporated, which we do here by using ZFITTERdH. 
In evaluating this expression to order we let can let both cw and sw take their on-shell 
values. 


The on-shell expression for which enters the Z couplings is that defined in Eq. ^1 
and can be written as 



m 


2 ( 0 ) 

w 


m 


2 ( 0 ) 


(31) 


In terms of the physical Z mass and the SM prediction for the mass, we hnd the following 
shift in the on-shell dehnition of ^ —>• Syy = ^ -|- 6 s^ where 






( 32 ) 
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At this point we next examine the interactions of the Z with the various fermion 
fields. The SM expressions for a fermion’s vector and axial-vector couplings to the Z are 
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f = 

V 


SGfMl 

71 





SGjMl 

72 



(33) 


These couplings are now modified in several ways: {i) the definition of Gf has been changed 
as described above, {ii) the couplings of the lighter SM-like Z are shifted via the modification 
in {Hi) the physical Z now has additional components from the heavy gauge fields. When 
these effects are combined the resulting shifts in and are given by 


Sgl 


SGfMl 

~7T~ 


(h 1 9^9^ 

V2 ^ 4 (^2 +^2)2 





/.2_ 


W 







{gtS Zh + g't^B^) - g'tQdB^ 


1 


SGfMl 

72 


1 9^9? 

V 4 (^2 +^2)2 




h 

2 


/ 


{gtSz^ + gt^Bh 


(34) 


For b quarks, there is an additional shift arising from vertex corrections containing the extra 
heavy isosinglet quark as well as those involving the top since its couplings have been modified 
as discussed above. Expressing this as a shift in the left and right-handed couplings, with 
9v = i9L + 5'Ij)/ 2 and g^ = {g'}^ — g'l)/2, and keeping only terms which are quadratically 
enhanced by the top quark or heavy isosinglet quark masses, we find 


Sgl = 


G< 


47r272 L . ^ , 


mrrijXiXj [gih {m, mj) - g^^h (m*, mj) 


dd-^w glS^h (mi) \ - mf [glh (m^, mt) 


-9^2 {mt,mt) 


9z {^w ~ dv] 


h {mt) + g'lh {mt) 


(35) 


where gz = and the functions Ik are defined in the Appendix. The Xi denote the 

mixing angles which arise when deriving the charged Goldstone coupling to the physical t 
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and t' states; we find Xt = cos (6**), x*/ = —sin {9t). We have explicitly snbtracted off the SM 
contribntion in Eq. EH This shift reqnires no renormalization, and is therefore a prediction 
of this model. We hnd that the right-handed conpling receives no contribntions which 
are qnadratically enhanced by the masses of the t or t'. 

Using the above expressions, we can derive the shifts in the precision EW observables. 
We examine the following approximately nncorrelated set of of observables: Mw, T;, 

Ni,, Rb, Ab, Rc, and Ac. Here E; is the total leptonic width of the Z, is the effective 

leptonic weak mixing angle on the Z pole, W is the ‘nnmber of nentrinos’ (dehned throngh 
the invisible width of the Z), Rb and Rc are respectively the ratios of bottom and charm 
qnark widths over the total hadronic width, and Ab and Ac are respectively the polarization 
asymmetries for b and c qnarks. We will also dicnss the deviation indnced in the NnTeV 
measurement of the on-shell weak mixing angle na. Once the various coupling shifts are 
determined, it is rather easy to obtain the contributions to these observables. We hnd: 


6Tf 



6 s 


2,lep 

W,eff 


6N, 

6 R, 

6A, 


alSgi - g[6g\ 


2>{2 


(9^5g^±9X 

V 9f + gf 


^had6Tq — T q6T had 


gf 


had 

-gf 


[gf + gf)‘ 


{gl^g, 



gX) . 


(36) 


The NuTeV experiment measures the following ratio of neutral and charged current cross 
sections for both neutrino and anti-neutrino beams scattering off an isoscalar target: 


R- 


^'nc ^nc 


1 

2 




(37) 


where in obtaining the second equality we have included the tree-level SM prediction for R_. 
In the present model, the charged current contribution, which now includes the exchange of 
both and W^, can be expressed entirely in terms oiGf. This is unchanged in this model 
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since G/ is taken as an inpnt parameter in the electroweak fit. The neutral current piece is 
modihed by both the shifted Zff couplings discussed above and the exchange of the new 
heavy gauge bosons and Z^. Given the masses and couplings of these gauge bosons, it is 
straightforward to compute the required diagrams and hnd the shift in R_. The masses of 
Bh and Zh are given in Eq. IHl writing the interaction between these states and the fermion 
helds as we hnd 

sP = 9th/2 , gP = -9tTi/2 , 

gP = 9, {Ti - 2Q') /2 , 92.B = -g,T’l2 . (38) 

The exact expression for is not very enlightening so we do not include it here. R- can be 
generally written in the form p{l/2 — in our analysis we will set p = 1 and interpret 

the deviation as a shift in the on-shell weak mixing angle 

From the discussion above it is clear that there are only four independent parameters 
that we need to consider in our numerical hts: 5, sin6't and the ratios of couplings gt/g and 
g't/g'■ All of the shifts in the electroweak observables as well as the predicted masses for the 
new gauge bosons can be written in terms of these few parameters. 


3 Results 

In this section we present our numerical analysis of the experimental constraints on the 
Littlest Higgs, including those from both direct searches and precision measurements. We 
attempt to identify throughout the features of our analysis that are independent of the 
fermion couplings to the hypercharge bosons. 

The six new TeV-scale states predicted in this model are the four new gauge bosons, 
Bh, Zh, and the isosinglet quark t!, and the isotriplet 0. Limits on the masses and 
couplings of such particles have been set by direct searches at high-energy colliders such as 
the Tevatron. The strongest constraints arise from the lack of observation for the production 
of Bh- This can be understood from examining Eq. IHl the factor of 1/V^ in its mass, which 
arises from the normalization of the hypercharge generator dehned in jB], and the fact that 
p' ~ .3 < p in the SM, show that it is predicted to be quite light relative to the scale /. We 
present the theoretical predictions for as a function of the coupling ratio x = g[/g in 
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Fig.E for several values of f] uib^ can become as small as 400 GeV for / as large as 2 TeV 
and x ~ 1. The 95% CL bounds resulting from direct searches at the Tevatron are shown in 
Fig. El As expected, the strongest constraints are for x ~ 1, where reaches its minimum 
value; in this region, the data bounds / ^ 3.5 TeV. The constraints weaken when x % 1. 
However, as we will see later, x ~ 1 is the parameter choice preferred by the precision EW 
data; the combination of the direct and indirect constraints will require / > 3.5 TeV for all 
X. We summarize the constraints on / arising from searches for Bh, and Zh production 
at the Tevatron Run I in Tabled as claimed above, the most signihcant limits are obtained 
from Bh production. We note that the constraints from and Zh production depend only 
on the coupling of the fermions to the SU(2)i gauge bosons, and not on their coupling to 
the hypercharge bosons. 


X 

fm 

/ 


1 

3.48 TeV 

.85 TeV 

.81 TeV 

1/2 

2.33 

.50 

.53 

2 

3.23 

.65 

.75 


Table 1: 95% CL lower bounds on / for specihc values of x arising from Tevatron limits on 
Bh, W^, and Zh production; for Bh, x = gjg, while for and Zh, x = gt/g- 


It is interesting to imagine how these constraints would evolve if the Tevatron Run II 
fails to hnd a signal for new gauge boson production in Drell-Yan collisions. What ranges of 
X would now be allowed as / is varied? The answer to this question can be found in Table El 
We see that the sensitivity to / would be extended by about 30% in comparison to Run I 
which would most likely exclude the possibility that / < 5 TeV. 


/ 

2 

2.5 

3 

3.5 

4 

4.5 

5 

5.5 

X 

< 0.24 

< 0.26 

< 0.30 

< 0.43 

< 0.51 

< 0.61 

0.74-1.68 

all 


Table 2: Estimated allowed ranges of x for different values of / in TeV following a null search 
for new gauge bosons at the Tevatron Run II assuming an integrated luminosity of 10 fb~^. 
For / = 5 TeV, the disallowed range is shown. 

We now analyze the constraints arising from EW precision data. The new parameters 
introduced in this model are the scale / (or the ratio S = v/f), the couplings gt and ht, and 
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Figure 1: The theoretical predictions for the mass of the gauge boson Bh as a function of 
X = g^/g , the ratio of its coupling relative to the corresponding SM coupling, for several 
values of the scale /. 

the top quark sector mixing angle Ot, which enters both 5p and Z-pole 6-quark observables 
through the loop corrections discussed in the previous section. Rather than perform a 
combined £t to all four quantities, we study a series of two-dimensional fits for various 
combinations of parameters; this is sufficient to illustrate the essential physics. In deriving 
our constraints we follow the analysis of uni; we calibrate our fits by first varying the Higgs 
boson mass until we hnd the value that reproduces the 95% CL upper bound on itih 
presented in pp. We then use this reference to determine whether a given set of model 
parameters provides a good £t to the EW data. We say that a set of parameters is disallowed 
if the resulting y^ exceeds this reference value. We will hrst consider the observables 
Ti, iVy, Rb, Ab, Rc, and Ac, and later will include as determined from NuTeV. 

The error correlation between these measurements is quite small P], and will be neglected 
in our analysis. We use ZFITTER to derive the SM predictions for these quantities. 

We hrst set gt = g and g[ = g , and vary both sin [9t) and 5. The results of this £t for 
both mu = 115 GeV and mn = 200 GeV are shown in Fig. |21 the shaded region is allowed 
by the EW £t at the 95% GL, while the remainder is excluded. Although mu = 200 GeV is 
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0.0 0.5 1.0 1.5 2.0 

X 


Figure 2: The expected number of Drell-Yan events from the production of the boson 
during Tevatron Run I (110 pb“^) as a function of a; = g^/g , the ratio of its coupling relative 
to the corresponding SM coupling, for several values of /. The horizontal line represents the 
95% CL bound from CDF CH for this mechanism; the parameter region above this line is 
excluded. 

allowed in this model, unlike in the SM where rriH < 193 GeV at the 95% CL, only a small 
sliver near the boundary (5 /a/ 2 = tan {9t) satisfies the EW precision constraints. This region 
persists but shrinks for larger values of mn- The tail along this boundary for large values 
of 6 is primarily due to the t and t' mixing contributions to 6p. We see that the Tevatron 
constraints are quite signihcant here, as they completely exclude this region. Included in 
these figures are curves representing the predictions for several t' masses as functions of 6 
and sin {6t) from Eq. 22. We see that the Tevatron limit requires that my > 7 TeV, which 
is significantly larger than the naturalness bound my < 2 TeV. This is beyond the range for 
direct discovery at the LHC. 

We now study the case where the couplings gt and g^ are varied within natural ranges 
away from their corresponding SM values. For this and all later discussions we assume 
ruH = 115 GeV. Shown in Fig. H] are the results of the EW fits for the following four 
coupling choices: (1) gt = g and g^ = g /2] (2) gt = g and g[ = 2g'] (3) gt = 5-/2 and g[ = g'] 
(4) gt = g/A and g't = g . We see from the top two hgures that the EW constraints are 
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0.20 


0.15 - 


0.10 


0.05 


0.00 



0.00 


0.05 0.10 

sin(0j 


0.15 



Figure 3: Fit to the EW precision data varying sin (6't) and 6 , for mn = 115 GeV (upper) 
and rriH = 200 GeV (lower). The diagonal line indicates the bound ( 5 / a /2 < tan(6't), the 
horizontal line denotes the 95% GL bound from production at the Tevatron. The series 
of curved lines indicates the t mass my as a function of 6 and sin {9t); from top to bottom, 
they represent my = 5, 7.5, 10, 15, and 20 TeV. The shaded regions are allowed by the EW 
£t. For the remainder of our discussion, we label g and g' as gsM and g'sM^ respectively, in 
the hgures. 
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strongest in the case where the Tevatron constraints are weakened. As can be seen from Fig. 
2, the Tevatron limits are the most stringent when g^jg k, 1, and decrease away from this 
region. This region is exactly where bpc vanishes, as can be seen from Eq. 1201 away from this 
point the contribution of to the EW £t becomes more important. For case (1), the EW 
precision constraints require / > 3.5 TeV, while in case (2) this bound is strengthened to 
/ > 7.2 TeV. These correspond to limits of > 7 TeV and > 14 TeV, respectively. The 
strongest constraints for cases (3) and (4) arise from the Tevatron direct searches; the EW 
precision constraints are relatively weak for these parameter choices. To further demonstrate 
that bounds from EW data become strong in precisely those regions in which the Tevatron 
constraints decrease, we present in Fig. Elthe results of hts where we set / = 2 TeV and vary 
both sin (6*^) and x. The three hts we perform use the following dehnitions of x: (1) gt = xg 
and g[ = xg'] (2) gt = xg/2 and g't = xg'] (3) gt = xg and g't = xg'/2. These variations in the 
couplings within natural ranges provide for a thorough exploration of the model parameter 
space. In all three cases the parameter values allowed by the EW fit are clustered very near 
g't = g', and are excluded by the Tevatron constraints. We obtain identical results if we 
choose / = 3 TeV. Combining the results of these hts, we can conclude that / > 3.5 TeV 
throughout the full parameter space of the model. This implies a minimum mass of my = 7 
TeV for the t' state, which is beyond the reach of the LHC. 

Having established that / ^ 3.5 is excluded by the data, we now set / = 4 TeV, which 
is allowed by both EW precision constraints and Tevatron direct search limits. We study the 
ht to the EW precision data as we vary both sin {6t) and x; we choose the same three cases 
for the relations of the couplings as examined above. The results of these hts are presented 
in Fig. ini We see that a far larger region of the coupling space is consistent with the data 
in this case, and that the EW data allows signihcant deviations from the g't = g limit. To 
examine the contributions of the various observables to the ht, we present in Fig. [7| the 
values for the various observables obtained in the gt = xg, g't = xg' ht; Ay^ is the diherence 
between the y^ values obtained in the Littlest Higgs model and in the SM. Included in these 
figures are all values of x which satisfy the bounds imposed by the EW data. We have not 
presented our results for the observables Rc, Ac, and because their Ay^ values are very 
nearly zero throughout the entire parameter space; i.e., they have no inhuence on the ht. 
All of the Ay^ ranges exhibit a dip either upwards or downwards at the right edge of these 
hgures. These arise from the t and t' contributions to either 5p or, for Rb, the Zbb vertex. 
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Figure 4: Fit to the EW precision data varying sin (9t) and 5, for various choices of gt and 
g^r as indicated. The diagonal line represents the bound 5/\/2 < tan (^t), the horizontal line 
denotes the 95% CL bound from Bh production at the Tevatron, and the series of curved 
lines corresponds to the t mass my as a function of 6 and sin (0^); from top to bottom, they 
represent my = 5, 7.5, 10, 15, and 20 TeV. The shaded regions are allowed by the EW £t. 
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m„ = 115 GeV 
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1 = 2 TeV 
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sin(et) 

Figure 5: Fit to the EW precision data varying both sin {6t) and x, for / = 2 TeV. The 
shaded regions are allowed by the EW £t, while the areas above the horizontal lines are 
forbidden by Tevatron direct search constraints. The rightmost boundary of the figure is 
determined by the constraint h/v^ < tan {6t). 


h EXCLUDED 
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gi=x*gsi„ g't=x*g'sj,/2 
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Tevatron 
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Although m^i becomes large as sin (6*t) —>■ 5l\/2, since the t! mixes with the top quark it does 
not decouple, in the sense discussed above, for some parameter space regions. We see that 
Mw is generally better predicted by the Littlest Higgs model than the SM; this remains true 
for other parameter choices. The shift in Mw depends only on the coupling of the fermions 
to the SU{2)i gauge bosons, and not on their coupling to the hypercharge bosons; this will 
therefore remain true for other choices of fermion transformation assignments. We also see 
that while Ri, is not affected until sin (0^) —> 5 j\/2, both T; and are very sensitive to 

the deviations predicted by the Littlest Higgs, and that the agreement between the predicted 
values for these observables and the experimental measurements is always worse than in the 
SM. 


Finally, we discuss the effect on the fit of the NuTeV measurement of the on-shell 
value of the weak mixing angle The NuTeV result currently disagrees with that derived 

from Z-pole data by approximately 3 standard deviations USI. Although this is possibly a 
signal of new physics, several more mundane explanations, such as larger parton distribution 
function uncertainties than those assumed by the NuTeV collaboration, asymmetric strange 
and charm sea quark distributions, or nuclear effects, have been suggested in|. With this 
caveat in mind, we will assume both the central value and error determined by the NuTeV 
collaboration, and repeat our analysis including the NuTeV result. To demonstrate that the 
allowed range of / is unaffected by the inclusion of the NuTeV measurement, we display the 
results of re-performing the following two analyses in Fig. |H| (1) varying 6 and sin {6t) while 
setting Qt = g and g[ = g' , the analog of the fit presented in Fig. 01 (2) varying sin {9t) and 
X with / = 2 TeV, and x defined by gt = xg and g[ = xg , the analog of the result shown in 
Fig. 01 Comparing these figures with those presented previously, we see that the constraint 
on / is unchanged by the addition of the NuTeV data. Although the full EW data set 
allows a slightly larger range of parameters here than before, the Tevatron constraint again 
excludes these regions. For other coupling choices, the NuTeV measurement again either 
very slightly increases or decreases the size of the allowed parameter space; no large swath 
of new parameter values are allowed. 

We next set / = 4 TeV, gt = xg and g't = xg', and vary both sin (6*^) and x; the 
results are presented in Fig. 01 We display the allowed region in the x — sin(6't) plane as well 
as the values of for The Littlest Higgs model shifts in the direction preferred 

by the NuTeV measurement throughout most of the entire parameter space; this remains 
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XI 



XI 



Figure 6: Fit to the EW precision data varying both sin {9t) and x, for f = A TeV. The shaded 
regions are allowed by the EW £t. The rightmost boundary of the figure is determined by 
the constraint (5/v^ < tan(0t). 
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Figure 7: The Ax^ values for several observables as functions of sin {6t) for the case gt = xg 
and g^ = xg ] included are all x values that satisfy the precision EW constraints as shown in 
the previous figure. 
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Figure 8: Fit to the EW precision data including the NuTeV measurement of the on-shell 
weak mixing angle. The top hgure studies the variation with 6 and sin (6*t), while the lower 
hgure shows the dependence on sin {9t) and x with / = 2 TeV. 
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true for other choices of the couplings. The size of the allowed region is smaller than in the 
case where the NuTeV data is not included. The shifts induced by the Littlest Higgs model 
(~ 0.5 — 1.5 a) cannot compensate for the large ~ 3 a deviation of the SM prediction from 
the NuTeV results. 

Before concluding we note that we have not discussed the possible influence of mea¬ 
surements of Atomic Parity Violation (APV) on the hts. Until recently the experimental 
value of the weak charge, Qw, differed from the predictions of the SM by about 2a PI- After 
a series of detailed theoretical calculations including various atomic effects it now seems that 
experiment and theoretical predictions are in agreement PI- However, since the theoretical 
calculation of Qw may still be in flux we do not include this observable in our hts. After 
computing the corrections to Qiv in this model and scanning the allowed regions from our 
previous hts, we hnd that either sign for the shift of Qiv can be generated. 


4 Conclusions 

Little Higgs models are a novel attempt to address the hierarchy problem. These theories 
predict the existence of a host of new particles at the TeV scale, which are necessary to cancel 
the quadratic divergences of the Standard Model. Here, we have examined the simplest 
version of these models, the Littlest Higgs, which contains four new gauge bosons, a weak 
isosinglet quark, t', with Q = 2/3, as well as an isotriplet scalar held with their masses being 
constrained by naturalness requirements. In this paper, we have considered the contributions 
of these new states to precision electroweak observables, and have examined the additional 
constraints provided by the direct searches for new particles at the Tevatron. We have 
performed a thorough exploration of the parameter space of this model and have found that 
there are small regions allowed by the precision data where the parameters take on their 
natural values. When the additional limits provided by the Tevaton are included, these 
regions are no longer permitted. By combining the direct and indirect effects of these new 
states we have constrained the ‘decay constant’/ to rather large values > 3.5 TeV; similarly, 
to satisfy both sets of data rrit' ^ 7 TeV which is far beyond the reach of direct searches 
at the LHC. These bounds imply that signihcant hne-tnning must be present in order for 
this model to resolve the hierarchy. We thus hnd that the Littlest Higgs model is tightly 
constrained by the combination of precision electroweak data and direct Tevatron searches. 
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NuTeV 



NuTeV 



sin(0j 

Figure 9: Fit to the EW precision data including the NuTeV measurement of the on-shell 
weak mixing angle; we have set / = 4 TeV and gt = xg, = xg . The top hgure shows 
the £t to the EW data, while the lower hgure illustrates the shift from the SM for those 
values of x and sin (dt) which satisfy the EW constraints. 
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Appendix A 


Here we investigate the effects of the scalar triplet vev on the precision electroweak 

observables. The vev gives a negative contribution to the p-parameter which can be written 
as 


V0 


^ r 

4 _a(p2 + p2)2 


(39) 


where A is the Higgs quartic coupling, which can be expressed in terms of the Higgs mass 

via [21] 


1 r ruh 

3 [200 GeV. 


(40) 


a represents a new parameter arising from the coefficient of an ultra-violet operator; natu¬ 
ralness suggests that the value of a lies in the range 0.1 — 1.0. Including this contribution 
in the fit to the global electroweak data set (minus the NuTeV measurement) yields the 
results displayed in Fig. (TUI Here, we have examined the previously least restrictive case, 
as observed from Fig. 4, and set the values of the couplings to gt = p/4, g[ = g', while 
varying the parameter a. For comparison, we also show our previous results for the case 
when the contributions from the triplet vev are not included. We see that the addition of 
the contribution from the triplet vev only worsens the ability of this model to accomodate 

the electroweak data set. We have checked that the constraints become even more restrictive 
as a approaches unity, and as the value of the Higgs mass increases. 


Appendix B 


We define here the loop integrals in Eq. |2H1 which arise in the computation of the shift 
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Figure 10: Results from the global fit to the electroweak precision data, including the con¬ 
tributions from the scalar triplet vev to the electroweak observables for various values of the 
parameter a. Also shown are the results from the fit when the triplet vev is not included. 
The horizontal line denotes the 95% CL bound from Bh production at the Tevatron, and the 
series of curved lines corresponds to the t' mass mt' as a function of 6 and sin(6'i); from top 
to bottom they represent = 2.5, 5, 7.5, and 10 TeV. The shaded regions are allowed by 
the EW fit. 
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in dg'l due to the vertex correction diagrams containing t and t' states: 


I2 


h ijrii) 


h {jrii) 


miiTij 


pi pl—x 

lo Jo {1 — X — y) — xyMz + xmj + yrn?^ 


pi pl—x 

dx dy\n 

Jo Jo 


{1 — X — y) — xyM^ + xrrif + ym 




pi pl—x 

I dx J dy\n 


(1 - y) -x{l-x-y)Ml + ymj 




nl 


dxxln 


xM^ + (1 — x) m. 




(41) 


where y is an arbitrary renormalization scale that cancels from the hnal expression after all 
the individual diagrams are summed over. 
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